KOI- 126: A Triply-Eclipsing Hierarchical Triple 
with Two Low-Mass Stars 
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The Kepler spacecraft has been monitoring the light from 150,000 stars in its 
primary quest to detect transiting exoplanets. Here we report on the detection 
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of an eclipsing stellar hierarchical triple, identified in the Kepler photometry. 
KOI-126 (A, (B, C)), is composed of a low-mass binary (masses M b = .2413 ± 
0.0030 Mq, Mc = 0.2127 ± 0.0026 M©; radii Rb = 0.2543 ± 0.0014 Rq, Rc = 
0.2318±0.0013 Rq-, orbital period Pi = 1.76713±0.00019 days) on an eccentric 
orbit about a third star (mass Ma = 1.347 ± 0.032 Mq-, radius Ra = 2.0254 ± 
0.0098 Rq; period of orbit around the low-mass binary Pa = 33.9214 ± 0.0013 
days; eccentricity of that orbit 62 = 0.3043 ± 0.0024). The low-mass pair 
probe the poorly sampled fuUy-convective stellar domain offering a crucial 
benchmark for theoretical stellar models. 

The Kepler mission's primary science goal is to determine the frequency of Earth-size plan- 
ets around Sun-like stars. To accomplish this, thousands of stars are monitored to detect the 
characteristic dimming of star light associated with a planet passing in front of its host (1). 

The Kepler observing scheme lends itself not just to planetary transits but, additionally, to 
the detection of the analogous eclipses in stellar multiple systems (2). We report on one such 
hierarchical system consisting of two closely-orbiting low-mass stars (KOT126 B, C; Pi ~ 
1.767 days) that are together orbiting with a longer period (P2 ~ 33.92 days) a more massive 
and much more luminous star (KOT126 A). The two low-mass stars double the number of fuUy- 
convective stars with mass and radius determinations better than a few percent; the stars in the 
eclipsing binary CM Draconis (3-5) are the remaining entries in this inventory. 

The KOI-126 system is oriented such that all three components and the Kepler spacecraft 
lie near a common plane and, consequently, eclipses among all three stars are observed. The 
unique variety of these eclipses reveal complex, information-laden dynamical content. As a 
consequence, the absolute system parameters for all three stars may be accurately determined 
from the photometry alone and are immune to the biases that plague the traditional study of 
low-mass eclipsing systems (e.g., as a result of stellar spots on the low-mass stars, 6). 
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The Kepler photometry, reported here, spans 418 days of near continuous observation of 
KOI-126 (KIC 5897826, 2MASS J 19495420+4 1065 14, Mag^ = 13.044) . The first 171 days 
consist of successive 29.426 min exposures, whereas the remaining 247 days were observed at 
both the 29.426 min (or long) cadence and at the 58.85 sec (or short) cadence (Fig. SI shows the 
long cadence time-series). We used the short cadence data (7) exclusively in our determination 
of system parameters. 

KOI-126 was originally identified by the Kepler data pipeline (8, 9) as a planetary candidate. 
The suspected transit events (Fig. 1) were found to be unusual upon closer inspection. The 1- 
3% decrements in the relative flux, occurring approximately every 34 days, exhibit a variable 
eclipse profile. At some epochs, two resolved transits are seen; at others, the transits are nearly 
simultaneous. The transit-like events have durations that vary substantially across the observed 
epochs and show asymmetric features about mid-eclipse, indicative of accelerations transverse 
to the line-of-sight. 

The periodic, superposed transits in the Kepler light curve are most readily explained as the 
passage of a close (or inner) binary (KOI-126 B, C) across the face of a mutually orbited star 
(KOI-126 A). At the start of each passage, KOI-126 B and C are at a unique phase in their 
binary orbit, yielding a unique transit route and light curve shape. To first order, the duration of 
transit is shorter or longer depending on whether the motion of a given component of the inner 
binary is prograde (shorter duration) or retrograde (longer duration) relative to the orbit of the 
inner binary center-of-mass. The short- timescale orbital motion of the inner binary accounts for 
the apparent accelerations. 

Preliminary modeling of the Kepler light curve predicted the secondary passage of the inner 
binary behind KOI- 126 A (shown in the bottom two time-series in Fig. 1). A Box-Least-Square 
algorithm (10) search of the Kepler light curve, excluding these secondary and transit events, 
revealed the relatively shallow eclipses between KOI-126 B and C, occurring every ^ 0.88 days 
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(Figs. S2, S3). Based upon these detections, KOI- 126 B and C were inferred to be each less 
luminous by a factor of ~3,000-5,000, as observed in the wide Kepler bandpass (1), than KOI- 
126 A. Both the eclipses of the inner binary pair and occultations of that pair by KOI- 126 A 
were not observed to be strictly periodic as they were absent in the data for long stretches of 
time. The long cadence event near t = 2, 455, 069.113 (BID) features the alignment of all three 
objects along the line-of-sight, resulting in a short brightening in the light curve (11). 

A periodogram of the light curve, after removing eclipse events and correcting for instru- 
mental systematics, shows an :^17 day modulation with a relative amplitude of ^500 parts-per- 
million (Fig. S4 plots a representative sample of this variation). It is likely associated with the 
rotation of KOI- 126 A (see Supporting Online Material, SOM). 

In addition to the Kepler photometry, we collected sixteen spectra of KOI- 126 over 500 
days (SOM). The spectra only showed features associated with a single star, KOI- 126 A. The 
primary goal of these observations was to acquire precision radial velocity (RV) measurements 
of KOI-126 A (Fig. S5, Table SI). An ancillary goal of the spectroscopic study was to fit model 
spectra to the composite spectra for KOI-126 A in order to determine stellar parameters (SOM, 
Table 1). From this analysis, we found that KOI- 126 A is metal rich, relative to the Sun ([Fe/H] 
= 0.15 ± 0.08), and has an effective temperature of 5875 ± 100 K and a stellar surface gravity 
log c/A = 3.94 ±0.14. 

We used a full dynamical-photometric model to explain the data (SOM). Newton's equations 
of motion along with a general relativistic correction to the orbit of the inner binary (12-14) 
were integrated to determine the positions and velocities of the bodies at a selected time. Each 
individual star's position was then corrected to account for the finite speed of light (77). These 
corrected positions, along with the absolute object radii and relative flux contributions, were 
used to calculate the combined flux. 

The model was fitted to the short cadence Kepler data by performing a least-squares mini- 
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mization (SOM). Only the data shown in Fig. 1 were utilized in the fit. While not used in the 
fit, the long cadence data are nearly exactly matched when using the best-fit parameters (Fig. 
S6). Subject only to the short cadence photometric data, we determine the individual masses 
and radii with fractional uncertainties less than 10% and 3%, respectively. 

We included the RV data for KOI- 1 26 A in a subsequent fit (Table 1 ). The best- fit parameters 
are identical to those found using photometry alone; although, masses were determined to better 
than 3% and radii to better than 0.5%. The inner binary orbit is nearly circular and inclined by 
PS 9° relative to the outer binary orbit. 

The photometric data could not be fit by assuming fixed Keplerian orbits for the inner and 
outer binaries. This is due to the relatively rapid variation of the Keplerian orbital elements 
(such as the orbital inclination or the eccentricity) as a result of the gravitational interaction 
between the three stars (Fig. S7) {12, 15). The observational evidence of this periodicity is 
most easily seen with the rapid circulation of the ascending node of the inclined inner binary: 
the inner binary orbit precesses, like a spinning top, every ~ 1000 days in response to the gravity 
of KOT126 A. This precession explains the occasional absence of the eclipses between KOI- 
126 B and C. Similarly, the oscillation of the outer binary inclination, compounded with the 
measured eccentricity, explains the sparsity of occultations of the inner binary by KOT126 A. 

Our final dynamical model did not include the effects of stellar spin, tidal distortion, or fric- 
tional dissipation because the data did not demand them. Nevertheless, we did investigate these 
contributions by extending our numerical model to include parameterized forces appropriate 
for each effect (SOM). The only relevant force on the observed timescale, assuming plausible 
rotation and dissipation rates, is caused by the mutual tidal distortion of KOT126 B and C. The 
controlling parameter for this force is the internal structure (or the apsidal) constant, k2, which 
is intimately related to the interior stellar density profile and provides an important constraint 
on stellar models {16, 17). There are no reliable constraints on k2 for stars with masses similar 
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to KOI- 126 B or C. For a fuUy-convective star, approximated as a polytrope with index n = 1.5, 
k2 is estimated to be as large as 0.15 (18). 

Our data constrain k2 < 0.6 at 95% confidence, assuming that KOl-126 B and C have 
equivalent apsidal constants. Although the current constraints on k2 are modest, a fit to the 
predicted Kepler light curve over the remainder of the nominal mission (3.5 years in total, 
with 25 more transit-like events) demonstrates that k2 ~ 0.1 will be measured with a relative 
precision of ~ 1%. In addition, masses and radii will be determined to better than 0.1%. 

The masses of KOI- 1 26 B and C lie below the threshold for having fully convective interiors 
(less than ^ O.SMq). It has been suggested (5, 19-22) that the large (10-15%) disagreements 
seen between model predictions (23, 24) and measured radii for low-mass stars were confined 
to stars in close binaries outside of this convective domain. A sample of low-mass stars with 
dynamically-measured properties validates this claim (3, 19-22, 25, 26); however, there is little 
reliable information available for stars with masses under O.SM© (see Fig. 2). 

Previously, CM Draconis (3-5) provided the only precise constraints on stellar models for 
stars below the convective mass boundary. In its case, theoretical models seemed to underesti- 
mate the stellar radii at the 5-7% level, a disagreement less than that seen with more massive 
stars but still consistent with increased activity attributed to fast rotation as a result of tidal spin- 
up (5). In comparison, the radii of KOI- 126 B and C are also underestimated by the models; 
however, this disagreement is smaller (2-5%; Fig. 2). 

In addition to accurately measured masses and radii, the metallicities and ages of KOI- 126 B 
and C are approximated by the values estimated for KOI- 126 A, if we assume all components 
were co-evolved and formed from the same proto-stellar nebula. In this case, KOI- 126 B and 
C have a super-solar metallicity ([Fe/H] = 0.15 ± 0.08) which can be compared to the poorly 
determined, sub-solar metallicity estimated for CM Dra (5). The enhanced metallicity of KOI- 
126 relative to solar may partially resolve the discrepancy between the observed and predicted 
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radii of KOI- 126 B and C, although we are not aware of any models that properly account for 
this enhancement. We compared the mass and radius of KOI- 126 A with stellar models in a 
well-calibrated mass range (27) and estimated a system age of 4 ± 1 Gyr. 

We are unable to measure the spin periods of KOI- 126 B and C, but, it is likely that synchro- 
nization has occurred and that the spin periods are nearly equivalent to the orbital period of the 
inner binary (SOM). These spin periods are slower than the expected orbit- synchronized spin 
periods for CM Dra A and B 1.3 days) by nearly 0.5 days. This fact may partially account 
for the differences in radii between the similar-mass stars CM Dra B and KOI- 126 C - CM Dra 
B may have increased magnetic activity relative to KOI- 126 C owing to its faster rotation (5). 
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Figure 1: Best-fit dynamical-pliotometric model to the fitted data for KOI- 126. The top 
eight panels are data observed during the passage of B and C in front of A from the perspective 
of the Kepler spacecraft. The bottom two panels are data observed during the passage of B 
and C behind A (binned by a factor of 10). In all cases, the solid black line gives the best- 
fit model. The inlaid diagrams show, to scale, the orbits of B (green) and C (blue) relative 
to A (yellow star) at times corresponding to those in the associated time-series. The dashed 
orbit is that of the center-of-mass of B and C. The numbered circles give the locations of B 
and C at the times indicated by the vertical dotted lines in the associated time-series, corre- 
sponding to the numeric index (0-2). The circles are to scale for the radii of B and C. The 
specific values To for each respective panel, reading from left to right and top to bottom, are 
(in BJD) 2455102.815, 2455136.716, 2455170.465, 2455204.267, 2455238.207, 2455271.751, 
2455305.713, 2455339.496, 2455259.000, 2455326.506. 

Figure 2: Masses and radii of known low-mass stars with dynamically-estimated masses 
and radii measured to better than 3% fractional accuracy (22, 28). The black curves correspond 
to the theoretical stellar isochrones by Baraffe et al. The dashed, dotted, solid and dash-dotted 
curves correspond to 1 Gyr, 2 Gyr, 4 Gyr and 5 Gyr solar metallicity isochrones, respectively. 
The blue points correspond to CM Dra A, B and the red points correspond to KOI- 126 B, C. 
The inset panel corresponds to the region in the larger plot enclosed by the dashed rectangle. 
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Figure 1 
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Figure 2 




Parameter 


Value 


Masses 




Ma 


1.347 ± 0.032 Mq 


Mb 


0.2413 ± 0.0030 Mq 


Mc 


0.2127 ± 0.0026 M0 


Radii 




Ra 


2.0254 ± 0.0098 Rq 


Rb 


0.2543 ± 0.0014 i?0 


Rc 


0.2318 ± 0.0013 i?0 


Average Densities (g cm~'^) 




PA 


0.2288 ± 0.0029 


Pb 


20.70 ±0.19 


PC 


24.09 ±0.20 


Surface gravities (logarithms in cgs units) 




log Qa 


3.9547 ±0.0069 


log gB 


5.0101 ± 0.0029 


log gc 


5.0358 ± 0.0027 


Observed relative fluxes 




Fa 


= 1 


Fb 


(3.26 ±0.24) X 10"'' 


Fc 


(2.24 ±0.48) X 10"" 


"Outer" binary [(A, (B, C))] elements on 2,455,170.5 (BJD) 




Period, P2 (day) 


33.9214 ±0.0013 


Semi-major Axis, a-z (AU) 


0.2495 ±0.0017 


Eccentricity, 62 


0.3043 ± 0.0024 


Argument of Periapse, uj2 


52.88° ±0.33 


Mean Anomaly, M2 


19.87° ±0.29 


Sky-plane Inclination, iz 


92.100° ±0.016 


Longitude of Ascending Node, fl2 


= 0° 


"Inner" binary [(B, C)] elements on 2,455,170.5 (BJD) 




Period, Pi (day) 


1.76713 ±0.00019 


Semi-major Axis, oi (AU) 


0.021986 ±0.000090 


Eccentricity, ei 


0.02234 ±0.00036 


Argument of Periapse, oji 


89.52° ± 0.42 


Mean Anomaly, Mi 


355.66° ± 0.42 


Sky-plane Inclination, ii 


96.907° ± 0.044 


Longitude of Ascending Node, fli 


8.012° ±0.039 


Star A parameters from spectroscopy 




Effective Temperature, Tcs^a (Kelvin) 


5875 ± 100 


Metallicity, [Fe/H] 


0.15 ±0.08 


Projected Rotational Velocity, va sini^ (km s"^) 


4.6 ±0.9 



Table 1: Parameters for KOT126 (A, (B, C)) based on a fit to the Kepler short cadence pho- 
tometry and sixteen radial velocity measurements. The final section of the table reports the 
parameters for KOI- 126 A based on an analysis of observed spectra having fixed the stellar 
gravity to that found from the fit to the photometry and RVs. 
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Online data All data used in this analysis have been made available at 

http : / / archive . stsci . edu/prepds/kepler_hlsp/ 

Spectroscopy Six of the sixteen spectra were collected using the TuU Coude Spectrograph on 
the 2.7m Harlan J. Smith telescope at the McDonald Observatory in west Texas, which has a 
resolving power of R=60,000 and a wavelength range of 3750-1 0000 angstroms . The remaining 
ten spectra were obtained using the Tillinghast Reflector Echelle Spectrograph (TRES; SI) on 
the 1.5m Tillinghast Reflector at the Fred L. Whipple Observatory on Mt. Hopkins, AZ. They 
were taken with the medium fiber, corresponding to R=44,000 and a wavelength range of 3850- 
9100 angstroms. 

The spectra were extracted and analyzed according to an established procedure (S2). We 
used multi-order cross-correlations to obtain precise relative velocities separately for the TRES 
and McDonald datasets. Spectral orders containing contaminating atmospheric lines were re- 
jected, along with orders in the blue with low signal-to-noise and orders in the red with known 
reduction problems. In total, we used 18 spectral orders in the McDonald cross-correlations, 
covering the wavelength range of 4500-6680 angstroms, and 1 8 spectral orders in the TRES 
cross-correlations, covering the wavelength range of 4580-6520 angstroms. The TRES veloc- 
ities were shifted to an absolute scale using the weighted mean offset from the single order 
velocities and the known TRES zeropoint (determined from the long term monintoring of lAU 
RV standards). The McDonald offset was fit for in the final joint solution. 

We shifted and coadded each dataset and classified the combined spectra by cross -correlating 
against the grid of Cf A synthetic spectra, which are based on Kurucz models calculated by John 
Laird and rely on a linelist compiled by Jon Morse. Interpolation between grid-points using 
the correlation peak heights yields the best- fit parameters, and the RMS scatter of the results 
from individual spectral orders provides an internal error estimate. Because of degeneracies in 
the stellar spectrum between T^s, loggA, and [Fe/H], these parameters are highly correlated, 
and an error in one parameter would likely result in systematic errors beyond the internal error 
estimates for the other parameters. For this reason, we double the errors to account for possible 
systematic effects. With log^^ fixed at 3.95, we find Teff=5875 ± 100 K, [Fe/H]=-i-0.15 ± 0.08, 
and va sin i=4.6 ± 0.9 km/s. 

We infer the rotational period of KOI-126 A to be 22±6 days, which is super-synchronous. 
We note that the observed 17 day periodicity in the Kepler light curve is comparable to the 
inferred rotational period, which suggests a possible association with the rotational modulation 
of surface features on KOI- 126 A. The variation is very roughly sinusoidal and out of phase 
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with the eclipse events. While the period of this variation is almost exactly half the period of 
the orbit of KOI- 126 (B, C) about KOI- 126 A, it is unlikely, given the amplitude and phase, that 
it is associated with the ellipsoidal distortion of A by the tidal field of (B, C) (S3). 

Dynamical-Photometric Model Positions and velocities. A hierarchical (or Jacobi) coordi- 
nate system is used when calculating the positions of the three bodies. In this system, ri is the 
position of C relative to B and r2 is the position of A relative to the center of mass of (B,C). We 
may specify ri and r2 in terms of osculating Keplerian orbital elements (period, eccentricity, 
argument of pericenter, inclination, longitude of the ascending node, and the mean anomaly: 
Pi,2, ei,2, ii,2, wi,2, ^^1,2, Mi_2, respectively). 

Newton's equations of motion, which depend on ri, r2 and the stellar masses, may be 
specified for the accelerations ri and r2 iS4, S5). An additional term may be added to the 
acceleration of ri due to the post-Newtonian potential of the inner binary (56). These are the 
only accelerations used in the fit giving the parameters listed in Table 1 . We worked in units 
such that G = 1. 

Further perturbing accelerations may be added to the acceleration of ri corresponding to 
the non-dissipative equilibrium tidal potential between B and C and the potential associated 
with the rotationally-induced oblate distortion of B and C (54). In this approximation, the axial 
spins of B and C follow the evolving orbit, staying normal to the orbit and spinning at a rate 
synchronous with the orbit. Both the accelerations due to tides and rotations depend on ri, Mb, 
Mc, Rb, Rc, k2,B and A;2,c- The acceleration due to rotation also depends on the angular axial 
spin rate of both B and C. The spin rates and apsidal constants are assumed to be the same for 
both B and C. We do not model the tidal or rotational distortion of A as their contributions to 
the total accelerations are negligible (S4). 

A final acceleration due to tidal damping may be added to the acceleration of ri (55). The 
scale of this acceleration is set by the tidal dissipation efficiency, Q. This acceleration is negli- 
gible for reasonable values of Q (> 100). 

We used an implementation of the Bulirsch-Stoer algorithm (57) to numerically integrate 
the coupled first-order differential equations for ri 2 and ri 2 in order to determine ri 2 and 
their temporal derivatives at any given time. The maximum step size in the integrator was 
Pi/1000 ^ 3 min. 

The Jacobi coordinates (ri and r2 and their derivatives) may be transformed into spatial 
coordinates of the three bodies, relative to barycenter. 

Radial velocity of KOI- 126 A. The RV data for KOI- 126 A were compared directly to the 
results of the numerical algorithm after applying a systematic offset associated with peculiar 
and bulk Galactic motion and an additional offset between the McDonald and TRES spectra to 
account for calibration error. The systematic offset was measured to be 7 = —27.278 ± 0.024 
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km s^^ and the additional offset was measured to be 7' = —0.26 ± 0.17 km s~^. 

Correcting for the finite speed of light. The positions of the three stars are projected to 
the location of the barycentric plane (i.e., the plane parallel to the sky-plane that includes the 
barycenter of KOI- 126 and that is normal to the line-of- sight) at a time to in order to correct for 
the delay resulting from the finite speed of light and the motion of the barycenter of KOI- 126 
along the line-of-sight. In detail, for each star and for a given observation (or "clock") time T, 
a secant line root finding algorithm is used to solve d{tr) — c{T — tr) — 7(T — to) = for the 
retarded time, t^., where d(t) is the star's distance at a time t from the barycentric plane, 7 is the 
radial velocity of the barycenter and c is the speed of light. The observed coordinates of each 
star at the clock time T are those found at each individual body's retarded time. For 7 = 0, 
the variation of the difference between the clock time and the retarded time is as much as a few 
minutes, in amplitude. We note that the light-time effects set the physical scale of the system 
and allow for the determination of absolute masses and radii {S8). Our code permits the use of 
a non-zero barycentric radial velocity 7 which may be measured from the radial velocity data. 
The effect of this parameter is to shift the model parameters with a fractional difference similar 
to 7/c ~ 10"^. The only parameters for which this effect is relevant are the orbital periods 
of the inner and outer binary which are adjusted by ^ 2.5 and ^ 0.88 sigma, respectively, 
relative to the periods measured with 7 = 0. All remaining parameters are adjusted by less 
than one tenth of their respective one sigma uncertainty. Given this fact and the possibility of a 
systematic bias in the velocity zeropoint, we opted to report (in Table 1) the parameters found 
assuming 7 = 0. 

Photometric model. The sky-plane projected 2D positions of all three objects were used as 
inputs to a light curve generating algorithm. 

All stars were assumed to be spherical. Additionally, the radial brightness profile of KOI- 
126 A was modeled as J(r)//(0) = 1 — ui(l — a/1 — r^) — ^2(1 — Vl — r'^Y where r is the 
projected distance from the center of A, normalized to the radius of A, and ui and U2 are the two 
quadratic limb-darkening parameters (59). The fluxes of B and C were specified relative to the 
flux of A. The sum of the fluxes was normalized to unity. The limb-darkening coefficients for 
A were found to be ui = 0.39 ± 0.03 and U2 = 0.22 ± 0.04. These values agree with theoretical 
expectations for stars with temperatures, gravities and metallicities similar to those found for 
KOI-126 A (0.3552 < ui < 0.4009 and 0.2553 < U2 < 0.2869 for 5750/i < Tes < 6000A', 
logfi'A = 4.0 and 0.10 <[M/H]< 0.20, SIO). We investigated various limb-darkened profiles 
for B and C according to the same model, however, any physically-valid choice of parameters 
gave similar fit values within the quoted uncertainties. 

For eclipses in which there are no three-way alignments, the loss of light is computed using 
an analytic prescription that depends on object separations, radii and relative fluxes (Sll). The 
total loss of light at a given time is the sum of the losses associated with the eclipse of A by B 



16 



or C, the eclipse of B or C by A, or the eclipses between B and C. When a three-way alignment 
occurs, the total loss of light is not given as a superposition of analytically-defined overlaps. In 
this case, the loss of light is computed numerically by integrating over the sky-plane. 

Determining best-fit parameters, covariances and uncertainties Fitting parameters. There 
were 23 free parameters in the final fit to all available data: the three masses (Ma, Mb, Mq), 
five orbital elements of the outer binary at to = 2,455, 170.5 (BID) {P2, 62, UJ2, ^2, M2), six 
orbital elements of the inner binary at to (Pi, ei, ui, ii,f2i,Mi), the radius of A (Ra), the 
relative radii of B and C (Rb/Ra, Rc/Ra), the relative fluxes of B and C (Fb/Fa, Fc/Fa), 
the two limb-darkening parameters of A {ui, U2), the systematic offset to the radial velocity 
of A (7) and a systematic offset between TRES and McDonald RVs (7'). The longitude of 
the ascending node of the outer binary is unconstrained and, for simplicity, has been fixed to 
= 0°. The longitude of the inner binary is measured relative to this orientation and does not 
reflect the true value. 

All short cadence data were initially fit with a multiplicative correction that was quadratic 
in time, to account for out-of-eclipse long-wavelength variability. The parameters describing 
these corrections correlated very weakly with the remaining 23 parameters and were therefore 
fixed to their best-fit values to reduce computation time. 

Best-fit parameters, errors and covariances. The data were fitted using an implementation 
of the Levenberg-Marquardt (L-M) algorithm (572, S13). The L-M algorithm minimizes the 
sum of the deviates squared (x^) by adjusting the model parameters. In addition to determin- 
ing the best-fit parameters, the L-M algorithm also returns the covariance matrix of the fitted 
parameters, approximating the x^-surface as quadratic in those parameters. For properly es- 
timated measurement uncertainty, the square root of the diagonal elements gives the formal 
1-sigma statistical errors in the parameters. The statistical error in any derived parameter may 
be estimated using the full covariance matrix. 

The measurement uncertainties for the photometric data were set equal to the root-mean- 
square deviation of the best-fit residuals for each short cadence event as plotted in Fig. 1 . The 
best-fit residuals were observed to show very little temporal correlation. The best-fit solution 
had a reduced-x^ of 1.00012 for 17370 degrees of freedom. The contribution from the RV 
data alone was 13.25. 

We note that while the masses Mb and Mc are known to no better than 1.5%, their ratio 
is known much more precisely: Mb/Mc = 1.1349 ± 0.0012. The ratio of radii are also well- 
determined: Rb/Ra = 0.12558 ± 0.00021 and Rc/Ra = 0.11445 ± 0.00015. Also, the 
sum of the radii of B and C are known relatively better than their difference - Rb + Rc = 
0.4861 ± 0.0026 Rq and Rb - Rc = 0.02254 ± 0.00026 Rq. The sum and differences of 
the masses of B and C have an analogous structure - Mb + Mc = 0.4540 ± 0.0057 Mq and 
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Mb- Mc = 0.02868 ± 0.00042 M©. 

We executed the L-M algorithm, starting from random positions in parameter space, many 
hundreds of times before finding the final solution. From this experience we note that a large 
number of local maxima (x^ minima) populate the likelihood landscape. In particular, solutions 
with inner binary periods belonging to a discrete set of aliases of the reported inner binary period 
satisfying P[ = P2/ {P2/P1 +n), for integers n, show qualitatively similar transit events. These 
solutions yield residual correlated structure, significantly higher and physically implausible 
stellar parameters for KOI-126 A. Moreover, these solutions do not predict the observed occul- 
tations or inner binary eclipses. A truncated Markov chain Monte Carlo search of the likelihood 
surface surrounding the adopted solution uncovered no larger extrema nor any local extrema 
with interestingly-low x^- 

Tidal Dissipation and Timescales Tidal dissipation is expected to modify the orbital ele- 
ments of the close binary and the spins of its components (514). In isolation from KOI-126 A, 
the orbital eccentricity of KOI-126 (B,C) should decrease with an exponential timescale ~ 1 
Gyr if Q/k2 = 10^ for each of B and C, where Q is the standard tidal quality factor (SI 5). In 
the presence of KOI- 126 A, the eccentricity would damp not to zero, but to a small fixed value, 
with the eccentricities of the close binary and third star aligned and precessing at the same 
rate (SI 6). However, the observed eccentricity ((ei) ^ 0.01) has not yet damped to that state 
(Figure S5). Moreover, if the system endured more than ~ 3 exponential damping timescales, 
then ei > 0.25 in the past, in violation of dynamical stability (SIT). Together with estimates of 
^2 — 0.1, these arguments bound Q > 10^. 

The nearly coplanar configuration of the triple system (9.2°, oscillating by 0.4° on each orbit 
of the third body) rules out significant eccentricity limit cycles (SI 8) or Kozai oscillations (SI 9) 
in the close binary. Aside from eccentricity damping, in isolation from KOI- 126 A, the spins 
of KOI- 126 B and C would damp quickly (SI 4) (~ 10*^ years) to an orbit-aligned and orbit- 
synchronized state. The presence of KOI-126 A causes the KOI-126 (B,C) orbit to precess 
faster than the natural spin precession of the component stars, so the spin vectors would damp 
to nearly the precession axis instead (S20). Although our observations are not currently sensitive 
to such spin-orbit misalignment, its long-term effect is to damp the mutual inclination (S19) on 
a timescale several to ten times longer than the eccentricity-damping timescale. 
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Figure SI: Long cadence Kepler light curve for KOI- 126 A showing data from Quarters 1-5. 
The plotted data have been normalized by a constant value but are otherwise the unprocessed 
product of Kepler aperture photometry. The discontinuities correspond to quarter breaks at 
which point the Kepler spacecraft is rolled and KOT126 is then observed at different pixels on 
the photometer. 
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Figure S2: Results of a search for eclipses between KOI- 126 B and C. Plotted is the "spec- 
trum" (dashed curve) from the Box-Least-Squares (BLS) search algorithm (S21) indicating pe- 
riodic dips in the Kepler photometry over 36 continuous days of observation, after excluding the 
prominent transit eclipses and de-trending by dividing by a 10 day moving average. The peak 
period corresponds to approximately half the orbital period of KOI- 126 B and C. The solid 
curve shows the result of the BLS algorithm after subtracting the best-fit photometric solution. 
A similar BLS search of data spanning an identical interval, observed 100 days later than this 
data, does not show any structure, indicating the eclipses had stopped. 
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Figure S3: Comparison of the observed and calculated light curve for the eclipses between 
KOI- 126 B and C. The observed data are those utilized in the BLS visually represented in Fig. 
S2. The calculated model is that using the best-fit parameters based upon a fit to the data in Fig. 
1 . The dashed line is 1:1. 
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Figure S4: Out-of-eclipse variability of KOI- 126. Plotted is a portion of the full Kepler light 
curve for KOI- 126 showing a ~ 17 day modulation. The solid vertical lines indicate the location 
in time of two transit events. 
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Figure S5: Radial velocity data for KOI- 126 A. The solid black curve gives the best-fit model 
to the data. The data with the larger (^0.5 km s~^) errors corresponds to McDonald measure- 
ments while the more precise data correspond to TRES measurements. 
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Figure S6: Long cadence data for KOI- 126 corresponding to the passage of KOI- 126 (B,C) in 
front of KOI- 126 A. The orange curve is the continuous model corresponding to the integrated 
best- fit model (in black). Refer to the caption of Fig. 1 for additional figure details. The specific 
values To for each respective panel, reading from left to right and top to bottom, are (in BID) 
2454967.760, 2455035.316, 2455069.113. 
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Figure S7: Evolution of the Keplerian orbital elements of KOI-126. The upper left panel shows 
the evolution of the eccentricity of the inner binary (composed of KOI- 1 26 B and C). The period 
of the secular variation in eccentricity is ~ 650 days. The large jumps in eccentricity correspond 
to periastron passage in the outer binary. The upper right panel shows the evolution of the inner 
binary inclination. The dashed lines bound the region in inclination for which eclipses of the 
inner binary will occur. The bottom left panel shows the evolution of the argument of periapse 
for the inner binary (dots) and for the outer binary [composed of the center of mass of KOT 
126 (B, C) and KOI-126 A; solid line] relative to the invariable plane that is perpendicular to 
the total angular momentum. The apsidal alignment of the inner and outer binaries corresponds 
to a maximum in the secular component of the inner binary eccentricity. The bottom right panel 
shows the evolution of the nodal longitude of the inner binary, relative to the invariable plane. 
The period secular variation of the inclination and the node is ~ 950 days. The span of time 
covered in the plots corresponds to the nominal 3.5 year Kepler mission. Times from 2,454,950 
(BJD) to 2,455,370 (BJD) cover the Kepler observations included in this work. 
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Time (BJD) Radial velocity of A (km s ^) RV Error (km s ^) 



2455339.898438 -26.651 0.101 

2455349.863281 -42.856 0.093 

2455366.925781 -4.787 0.087 

2455367.765625 -3.123 0.097 

2455456.621094 -34.408 0.061 

2455463.695312 -17.632 0.062 

2455466.714844 -9.171 0.063 

2455483.679688 -44.097 0.064 

2455485.667969 -41.976 0.073 

2455488.613281 -37.594 0.063 

2455051.753906 -34.411 0.400 

2455076.828125 -44.535 0.401 

2455311.859375 -44.997 0.401 

2455341.812500 -37.164 0.401 

2455394.703125 -21.844 0.401 

2455396.714844 -16.319 0.401 



Table S 1 : Radial velocity measurements for KOT 1 26 A. The adopted error is the sum in 
ture of the internal error and the estimated error attributed to instrument stability. 
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